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ABSTRACT 

If dark matter has a large self-interaction scattering cross section, then interactions among 
dark-matter particles will drive galaxy and cluster halos to become spherical in their centers. 
Work in the past has used this effect to rule out velocity-independent, elastic cross sections 
larger than a/m ~ 0.02 cm 2 /g based on comparisons to the shapes of galaxy cluster lens- 
ing potentials and X-ray isophotes. In this paper, we use cosmological simulations to show 
that these constraints were off by more than an order of magnitude because (a) they did not 
properly account for the fact that the observed ellipticity gets contributions from the triaxial 
mass distribution outside the core set by scatterings, (b) the scatter in axis ratios is large and 
(c) the core region retains more of its triaxial nature than estimated before. Including these 
effects properly shows that the same observations now allow dark matter self-interaction cross 
sections at least as large as a/m = 0.1 cm 2 /g. We show that constraints on self-interacting 
dark matter from strong-lensing clusters are likely to improve significantly in the near future, 
but possibly more via central densities and core sizes than halo shapes. 
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1 INTRODUCTION 

The nature of dark matter is one of the most compelling myster- 
ies of our time. On large scales, the behavior of dark matter is 
consistent with what cosmologists of yore called "dust" (e.g., |Tol^| 
|man|[T9 34), meaning its behavior is consistent with being colli- 
sionless and non-relativistic ("cold") for the vast majority of the 
Universe's history ( |Reid et al.|2010) . This consistency has been of 
great interest to the particle-physics community because the most 
popular candidate for dark matter, the supersymmetric neutralino, 
displays exactly this behavior iSteigman & Turner 1985, Griest| 



[T9i8l|Jungman, Kamionkowski & Griest|1996) . While the super- 
symmetric neutralino paradigm is attractive in many ways, there 
are two outstanding problems with it. First, astroparticle searches 
have yet to turn up evidence for the existence of the neutralino, 
though searches are rapidly increasing their sensitivity to interest- 
ing neutralino parameter space (Geringer-Sameth & Koushiappas 
20TT] [Ackermann et al.|[25TT] |Cotta et al.|2012HFox et al.|20TTT 
Bertone et al.|2012||Atlas Collaboration|2012||Koay|2012||Baudis| 

pra-| 



2012; 



Baer, Barger & Mustafayev||20 1 2 



IXENON100 Collabora- 



tion et al.|2012| l. Second, there are predictions for the structure of 



dark-matter halos that have not been observationally verified at a 
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quantitative level (Dubinski & Caiiber g|T99l]|Diemand et al.|2008| 
|Vogelsberger et al'|2009||Stadel et al.|2009||Navarro et al.|2010) . In 
fact, there are hints of tension with the neutralino paradigm on sub- 



galactic scales (Dobler & Keeton 2006; Gentile et al. 2007 


de Blok 


etal.|2008; 


Kuzio de Naray, McGaugh & de Blok 2008 


Cuzio de 


Naray et al. 


2010, Zwaan, Meyer & Staveley-Smith 2010; Boylan- 


Kolchin, Bullock & Kaplinghat 201 1 ; Papastergis et al. 201 1 ). And 



yet, most of the effort to characterize the evolution of the Universe 
(experimentally, observationally, and theoretically) has been in the 
context of this dust-like cold dark matter (CDM). In the absence of 
evidence for any dark-matter candidate, much less the neutralino, it 
is important to explore the structure and evolution of the Universe 
for dark-matter phenomenology beyond cold and collisionless. 

One intriguing possibility is that the dark matter belongs to 



and interacts with a "dark" or "hidden" sector ( jKhlopov, Stephan& 
Fargion|2006"l|Feldman, Kors & Nath|2007||Fo"ot|2007[|Pos pelov, 
Ritz & Voloshin|2008[|Feng & Kumar|2008[ |Arkani-Hamed et al. 
2009[|Feng et al.|2009||Sigurdson|2009| |Cohen et al.|2010J . the 



Standard Model has rich intra-sector phenomenology; it is not un- 
reasonable to speculate that a complex dark sector only tenuously 
connected to the Standard Model might exist as well. The sim- 
plest kind of dark-matter interaction is a hard-sphere interaction 
of identical dark-matter particles. Such an interaction - with an 
isotropic, velocity-independent, elastic scattering cross section - 
was first introduced in an astrophysical context by Spergel & Stein- 
|hardt| {2000[ >. Interactions of this type were invoked to ameliorate 
the tensions between observations and CDM predictions on small 
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scales (on the scales of individual dark-matter halos) while leaving 
the large-scale successes of CDM intact. In this paper we revisit this 
basic class of self-interacting dark matter (SIDM) using cosmolog- 
ical simulations to explore its effect on dark matter halo shapes as a 
function of cross section. In a companion paper (Rocha et al. 2012) 
we investigate implications for dark matter halo substructure and 
density profiles. 

We are reinvestigating this simple SIDM model, which had 
been decreed "uninteresting" in several studies a decade ago, for 
two primary reasons. First, we suspected that the constraints that 
indicated that the SIDM cross section was too small to meaning- 
fully alter the morphology of dark-matter halos, were not as tight 
as claimed. Second, there is a wealth of new data (e.g., from near- 
field cosmology, lensing studies of galaxies and clusters) that may 
be better places to either look for SIDM or constrain its proper- 
ties. In this paper and our companion paper, Rocha et al. (2012), 
we reevaluate past constraints on SIDM and suggest several new 
places to look for the effects of SIDM on halo structure. 

There was a burst of work on SIDM before its untimely demise 
I Yoshida et al.|2000a|b||Kochanek & White|2000| |Hogan & Dal 
canton||2000| IDalcanton & Hogan||2001||Dave et al.||2001| |Coirn 



Table 1: Summary of simulations. 



et al. 2002, Hennawi & Ostriker 2002). The death of isotropic, 
velocity-independent, elastically scattering SIDM largely came 
from the interpretation of three types of observations: halo evapo- 
ration in galaxy clusters ( Gnedin & Ostriker 2001 1, cores in galaxy 

and halo 



|200l} , 



clusters \ Yoshida et a l. 2000b; Meneghetti et al 
shapes (Miralda-Escude 200211. The Y2K-era constraints from the 
former two classes of observations were at the level of a/m < 
0.3 cm 2 /g and a/m < 0.1 cm 2 /g respectively. However, we show 
in a companion paper, Rocha et al. (2012), that the evaporation and 
cluster-core constraints are likely overestimated. 

The most stringent constraints on dark matter models with 
large isotropic, elastic self-scattering cross sections emerged from 
the shapes of dark-matter halos, in particular from lens modeling 
of the galaxy cluster MS 2137-23 by |Miralda-Escude| p002). This 



massive galaxy cluster has a number of radial and tangential arcs 
within ~ 200 kpc of t he halo center {M ellier, F ort & Kneib|1993| 
|Miralda-Escude|1995| >. |Miralda-Escude| l |2002} argued that self in- 
teractions should make dark-matter halos round within the radius r 
where the local per-particle scattering rate equals the Hubble rate 
r(r) = Ho, or equivalently, where each dark-matter particle ex- 
periences one interaction per Hubble time. The scattering rate per 
particle as a function of r in a halo scales in proportion to the local 
density and velocity dispersion, 



r(r) ~ p(r){o I 'm)v lms (r) 



(1) 



where p is the local dark-matter mass density and v Ilns is the rms 
speed of dark-matter particles. Using the fact that the lens model 
needs to be elliptical at 70 kpc, Miralda-Escude (2002) set a con- 
straint of a/m < 0.02 cm 2 /g on the velocity-independent elas- 
tic scattering cross section. This constraint is one to two orders 
of magnitude tighter than other typical constraints on velocity- 
indep endent scattering fYoshida et al.||2000bl |Gnedin & Ostriker] 
|2001| |Randall et al.|2008| >. It rendered velocity-independent scat- 
tering far too small to form cores in low surface brightness galaxies 
and other small gala xies |Kuzio de Naray, McGaugh & de Blok| 
2008; de Blok et al. 2008). This is unfortunate because the main 
reason SIDM was interesting at the time was that it was a mecha- 
nism to create cores in such galaxies (Spergel & Steinhardt 2000) . 

The tightness of the SIDM constraints on cluster scales has 
meant that the focus of SIDM studies has shifted to those on 
velocity-dependent cross sections, such that SIDM may signifi- 
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^Box 


m p 
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a/m 




[Mpc/h] 


[M /h] 


[kpc/h] 


[cm 2 /g] 


CDM 


50 


6.9 x 10 7 


1.0 






25 


8.6 x 10 6 


0,1 




SIDMq.i 


50 


6.9 x 10 7 


1.0 


0.1 




25 


8.6 x 10 6 


0,1 


0.1 


SIDMi 


50 


6.9 x 10 7 


1.0 


1.0 




25 


8.6 x 10 6 


0,1 


1.0 



*Note: columns give name, simulation box size, particle mass, force 
resolution, and interaction cross section. We use h = 0.71. See Rocha et 
al. (2012) for more details on the simulations. 

cantly alter dwarf-scale or smaller dark-matter halos while leaving 
cluster-mass halos largely untouched l Feng et al. 2009; Buckley 



& Fox|2010[|Loeb & Weiner|201 l| |Vogelsberger, Zavala & Loeb 
2012| >. In recent times, such velocity-dependent interactions have 
arisen in hidden-sector models designed to interpret some charged- 
particle cosmic-ray observations as evidence for dark-matter anni- 
hilation ([Pospelov, Ritz & Voloshin||2008l [FoT & Poppitz] [2009 



|Arkani-Hamed et al.||2009| |Feng, Kaplinghat & Yu||2010) . Con 



straints on other hidden-sector dark-matter models have been made 
using X-ray isophotes of the gas in the halo of the elliptical galaxy 



NGC 720 (Buote et al. 2002, Feng et al. 2009; Feng, Kaplinghat & 


Yu|2010||Buckley & Fox 


2010 Ibe&Yu|2010 McDermott, Yu & 


Zurek|2011 Feng, Renta 


a& Surujon|2012i. 



However, as we show below, reports of the death of isotropic, 
velocity-independent elastic SIDM are greatly exaggerated. In this 
paper, we show that these earlier studies did not correctly account 
for the fact that the observed ellipticity (of the mass in cylinders or 
the projected gravitational potential) gets contributions from mass 
well outside the core, and the region outside the core retains its 
triaxiality. We also show that for ellipticity estimators that are rel- 
evant observationally, there is significant amount of scatter and the 
overlap between CDM and SIDM ellipticities is substantial even for 
a/m = 1 cm 2 /g. Lastly, we find that in the regions where SIDM 
particles have suffered (on average) about one or more interactions, 
the residual triaxiality is larger than what has been previously esti- 
mated ( |Dave et al.|2001[ >. Along with the analysis in Rocha et al. 
(2012), we find that studies of the central densities of dark-matter 
halos are likely to yield tighter constraints on the SIDM cross sec- 
tion than the morphology of the halos. 

We briefly summarize our simulations in Sec. [2] We present 
results on the three-dimensional shapes of SIDM dark-matter ha- 
los compared to their CDM counterparts in Sec. [3] We reexam- 
ine the previous SIDM constraints based on halo shapes in light of 
our simulations in Sec. [4] In particular, we reexamine the Miralda-| 
|Escude| J2002) constraint in Sec. 4.2 and from the shapes of the 



X-ray isophotes of NGC 720 (Buote e t al.|2002{ l in Sec|43] In Sec. 
|4.4| we show how other lensing data sets may constrain SIDM in 
the future. We summarize the key points of this paper and present 
a few final thoughts in conclusion in Sec. [5] 



2 SIMULATIONS 

We modeled self interactions by direct-simulation Monte Carlo 
with a scattering algorithm derived in Appendix A of Rocha 
et al. (2012) and implemented within the GADGET-2 ( |Springel| 



SIDM Halo Shapes 3 



Major axis, CDM SIDM l SIDM X 




kpc/ h kpc I h kpc/h 



Inner region 



100 



-100 




100 



5, 

4i 



-100 




100 



*5i 



-100 




Figure 1. Surface density of a halo of mass M v j r = 1.2 X 1O 14 M0 projected along the major axis of the moment-of-inertia tensor - the orientation that 
dominates the lensing probability. The left column shows the halo for CDM, while the middle and right columns show the same halo simulated using SIDM 
with a/m = 0.1 cm 2 /g and 1.0 cm 2 /g, respectively. The bottom row shows the same information, now zoomed in on the central region. The surface density 
stretches logarithmically from R2 10~ 3 g/cm 2 (blue) to 10 g/cm 2 (red). 



|2005[ > cosmological N-body code. Once the code passed accuracy 
tests, we performed cosmological simulations of CDM and SIDM 
with identical initial conditions for cubic boxes of 25/i _1 Mpc 
on a side and 50ft _1 Mpc on a side, each with 512 3 parti- 
cles. For the SIDM runs we explored cross sections of a/m — 
1, 0.1, and 0.03 cm 2 /g, though the lowest cross section run (with 
a/m = 0.03 cm 2 /g) provided results that were so similar to CDM 
that we have not included them in any of the figures below. The ini- 
tial conditions were generated using the MUSIC code at z = 250 
( Hahn & Abel|2011) with a year seven WMAP cosmology ( |Ko-| 
|matsu et al.||2011) : h = 0.71, Q m = 0.266, fi A = 0.734, 
fib = 0.0449, n s = 0.963, a 8 = 0.801. A summary of our sim- 
ulation parameters, including particle mass and resolution is pro- 
vided in Table 1. We adopt a naming convention where the simu- 
lations denoted SIDMo.i and SIDMi have subscripts correspond- 
ing to their cross sections in units of cm 2 /g. In all cases the self- 
interaction smoothing length as defined in Rocha et al. (2012) was 
set to 2.8 times the force softening. 

We locate and characterize halos using the publicly available 
Amiga Halo Finder (AHF; Knollmann & Knebe 2009 ) package. 
The total mass of a host halo M v i r is determined as the mass within 
a radius r v k using the virial overdensity as defined in Bryan & 
|Norman| ( |1998| >. Though most of our analysis focuses on distinct 
(field) halos, we also explore subhalo shapes. For these objects their 
masses are measured within the radius at which the radial density 



profile of the subhalo begins to rise again because of the presence 
of the host. 



We present results on the shapes of halos at z = in a radial 
range r m i n to r v ir, where r m i n is the minimum radius within which 
we trust the shape measurements. Since our two sets of simula- 
tions have different resolution, we can check the convergence of our 
shape estimates. For integral measures of shape (e.g., the moment- 
of-inertia tensor for all particles within a given radius), we find that 
although the density profiles look largely converged outside of the 
numerical-relaxation radius r rc i ax defined in |Power et al.| (20031, 
the shapes do not converge until at least r m i n = 2r re iax- This ra- 
dius is roughly r m i n w 20 kpc for the halos in the 50/i _1 Mpc- 
sized simulations, with only modest dependence on halo mass and 
scattering cross section, and r m i n « 10 kpc for the 25/i _1 Mpc 
boxes. Below r m i n we find that the halo shapes are systematically 
too round. However, shapes are more robust if found in shells (ei- 
ther spherical or ellipsoidal) because they are less contaminated by 
the effects of the overly round and numerically relaxed inner re- 
gions. Shape estimates, especially integral estimates, are most re- 
liable if there are at least ~ 10 4 particles within the virial radius 
(or tidal radius for subhalos), consistent to what has been found in 
earlier work (e.g., |Allgood et al.|2006|[Vera-Ciro et al.|2011| >. 
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Figure 2. Surface density profiles for the same halo shown in Fig.[T] now projected along the intermediate axis. Deviations from axisymmetry are highest 
along this projection. 






Figure 3. Host halo shapes in shells of radius scaled by the virial radius in three virial-mass bins as indicated. The black solid lines denote the 20th percentile 
(lowest), median (middle), and 80th percentile (highest) value of c/a at fixed r / r v ; r for CDM. The blue dashed lines show the median and 20th/80th percentile 
ranges for <j/m = 1 cm 2 /g, and the green dotted lines show the same for cr/m = 0.1 cm 2 /g. There are 440, 65, and 50 halos in each mass bin (lowest mass 
bin to highest). 



3 SIMULATED HALO SHAPES 
3.1 Preliminary Illustration 

Before presenting a statistical comparison of CDM and SIDM halo 
populations, we provide a pictorial illustration of how an individ- 
ual halo changes shape as we vary the cross section. The columns 



of Figs. [T] and [2] show surface density maps for the same halo sim- 
ulated in CDM, SIDMo.i, and SIDMi from left to right. In Fig. 
[T] we project the halo along the major axis, which is the orienta- 
tion that maximizes the strong -lensing cross section (van de Ven, 
Mandelbaum & Keeton 2009, Mandelbaum, van de Ven & Kee- 
|ton||2009") . In Fig. [2] we project the halo along the intermediate 
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axis, which maximizes the deviation of the surface density from 
axisymmetry. This particular halo is one of the most massive halos 
identified in the 50/i _1 Mpc box runs, withMvir = 1.2 x 1O 14 M 
and r v i r = 1.27 Mpc. The top row in each figure shows the surface 
density on the scale of the virial radius, while the lower row shows 
the inner 300/i _1 kpc of the halo (side-to-side). The surface-density 
stretch is the same between the two figures. 

The major and minor axes for the projections in these figures 
were determined using the moment-of-inertia tensor of all parti- 
cles within a sphere of radius r v - 1T in the halo. If modeling the 
mass distribution as an ellipsoid, the principal axes a(major) > 
^(intermediate) > c(minor) are the square roots of the eigenval- 
ues of this tensor. 

In comparing Figs.[T|and[2] note that Fig.[T]is the most relevant 
for strong-lensing studies (Sec. |4j and shows the smallest differ- 
ences, especially at large radii. Indeed, only the zoomed view of the 
a/m = 1 cm 2 /g run is visibly rounder than the CDM case. Even 
for the intermediate projection (Fig.[2j, which maximizes the visual 
difference, the inner regions of the halo are only slightly rounder 
and less dense than their CDM counterparts for a/m — 0.1 cm 2 /g. 
The a/m — 1 cm 2 /g case is indeed less dense and rounder within 
~ 100/i -1 kpc, but even in this case, some ellipticity is clearly ev- 
ident. 

A final point of interest in these visualizations concerns the 
substructure. The subhalos apparent in the CDM halo are similarly 
abundant in the SIDM cases, and even approximately match in their 
positions. There are minor differences in substructure densities and 
locations (especially in the central regions) but overall it is difficult 
to distinguish among the runs by comparing their substructure con- 
tent (see Rocha et al. 2012 for a more quantitative comparison of 
substructure). 



3.2 Three-dimensional halo shapes 

We quantify halo shapes by examining ellipsoidal shells centered 
on the radial slices identified by AHF ( Knollmann & Knebe [2009} 
for profile measurements. We use shells instead of enclosed volume 
because it is less sensitive to numerical relaxation effects at the 
center and because it is a better estimate of the effects of local dark- 
matter scattering. In each shell of material, we calculate a modified 
moment-of-inertia tensor (defined and used in Allgood et al ,|2006} 
in ellipsoidal shells, 

shell 

ifj(a) = ^ ' ,T V'" » where (2) 

n 

r n = y/xl n +xlJ{b/ay+xlJ{c/ay 

and (a:i,n, X2, n , Xz, n ) are the coordinates of the nth particle in the 
frame of the principal axes (major, intermediate, minor) of this ten- 
sor. This is the same moment-of-inertia tensor from which shapes 
are inferred in |Dubinski & Carlberg|(l99l] > and |Dave et allpOOT} . 
The principal axes (a, b, c) are computed as the square roots of the 
eigenvalues of Xfj'.The weighting of the moment-of-inertia tensor 
is chosen such that the outermost particles in the shell do not domi- 
nate the shape estimate. We begin by finding the moment-of-inertia 
tensor in a spherical shell, setting a = b = c = 1 for this initial 
estimate of Xfj , and iterate to find Xff with convergent (a, b, c) 
values. In each iteration, the ellipsoidal shell volume is defined us- 
ing the (a, 6, c) found in the previous iteration. We experiment with 
either keeping the semi-major axis a of the shell fixed between it- 
erations or allowing a to float such that the volume in the shell 
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Figure 4. Median subhalo shape vs. radius for galaxy-mass systems com- 
pared to host halos of the same mass. Bold lines denote subhalos and lighter 
lines denote host halos. The radii are normalized by r v ir for hosts and T" t idai 
for subhalos. Line colors and styles have the same meanings as in Fig. [3] 
dashed blue is SIDMi, dotted green is SIDMo.i, and solid black is CDM. 

remains fixed as we iterate to find (a), but find that c/a is insen- 
sitive to these choices. Throughout this section, we show c/a for 
fixed volume in the shell. We only show results for c/a because the 
trends for b/a are similar but less informative, since c/a indicates 
the deviation of the halo shape from sphericity. 

In order to understand trends, we split our analysis of host 
dark-matter halos and subhalos, and bin halos by virial (or tidal) 
mass. Host halos are those whose centers do not lie within the 
virial radius of a more massive halo. In Fig. [3] we show the minor- 
to-major axis ratio c/a as a function of radius normalized by the 
halo virial radius, r/r v i r , for host halos in three mass bins. Larger 
values of c/a imply more spherical halos. For the two lower mass 
bins, we used halos selected from the 25 /i _1 Mpc boxes since these 
have the higher resolution. For the highest mass bin, we used halos 
identified in the 50 h~ 1 Mpc boxes in order to gain better statistics. 
We checked to make sure that the results were convergent between 
boxes where the relevant mass resolutions overlap. The shaded re- 
gion corresponds to the 20th to 80th percentile for c/a of the halo 
population for fixed r/r v i r , and the central line shows the median 
value of c/a. The black solid lines and yellow shaded regions de- 
note shapes of CDM halos, the blue dashed lines and regions cor- 
respond to a/m — 1 cm 2 /g, and the green dotted lines and re- 
gions correspond to a/m = 0.1 cm 2 /g. The regions extend down 
in r/rvir to the largest value of r m i n /r v i r in the given mass bin. 

We reproduce the well-known trend that galaxy-mass halos in 
CDM are more spherical than cluster mass halos and that CDM 
halos become more spherical in their outer parts ( Allgo od et al.| 
|2006[ >. SIDM halos deviate most strongly from CDM at smaller 
radii, where the scattering rates are highest for a fixed cross section. 
For SIDMi, halos are actually more spherical in their centers than 
their edges, with c/a rising with decreasing r for r / r v - 1T < 0.5. For 
SIDMo.i, differences from CDM are only apparent for r/rvir < 
0.1. 

In Fig. [4] we compare the shapes of subhalos (thick lines) and 
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host halos (thin lines) of similar mass by plotting the median axis 
ratio c/oasa function of r/r outcr where r out er = mdai for subha- 
los (as defined by the AHF halo finder) and r ou t cr = r v i r for host 
halos. All halos have masses within r out or between 10 11 h -1 Mq 
and 1O 12 /i _1 M . Though not shown, the 20th to 80th percentile 
ranges are similar in size as in Fig. [3] We find that the subhalo inte- 
riors in the a/m = 1 cm 2 /g cosmology are systematically rounder 
than host halos. We speculate that there are are at least three effects 
that drive this trend. First, subhalos are typically more evolved than 
field halos of the same mass, with fewer recent mergers. Allgood 
|et al.| ( |2006| l find that halos that form earlier are more spherical 
than halos that form later, which is attributed to directional merg- 
ing, and more generally to the highly non-spherically-symmetrical 
way in which halos form and accrete. Second, these subhalos are 
the remains of more massive halos, which are more susceptible to 
the effects of self interactions at fixed r/r v i r , as we showed in Fig. 
[3] Moreover, the outer radius of the subhalos is truncated with re- 
spect to its virial value, thus increasing r/r outcI for fixed r. Since 
halos tend to be rounder in the outer parts with respect to the in- 
teriors for CDM halos, this boosts the initial c/a for fixed r/r v i r , 
beyond which SIDM boosts c/a even more. We see this trend for 
the CDM and a/m = 0.1 cm 2 /g halos in Fig.|4] 

To see how the shape of dark-matter halos changes as a func- 
tion of the typical local scattering rate (Eq. [TJ, weplot c/a as a 
function of p(r)v TlnB (r) oc r(r)(a/m)~ 1 in Fig.pl The propor- 
tionality constant in relating pv lnla to the scattering rate is 0(1) 
and depends on the distribution function of dark matter particles. 
Thus it is reasonable to use pv ims as a proxy for the local scatter- 
ing rate modulo the actual cross section. To simplify the interpre- 
tation further, we multiply this quantity by 10 Gyrcm 2 /g in Fig. 
[5] In these units, if p(r)v lmB (r) > 1 for a/m = lcm 2 /g, most 
particles will have scattered after 10 Gyr. For a/m = 0.1cm 2 /g, 
this quantity needs to be 10 times larger to achieve the same scat- 
tering rate. Generally, pw rm s increases as one goes in towards the 
halo center, so particles tend to scatter more frequently in the core 
than in the outer parts of the halo, where interactions are uncom- 
mon over a Hubble time. 

Fig.[5]shows that deviations in the halo shape from CDM begin 
when pw rms (cr/m) x (10 Gyr) ~ 0.1, independent of halo mass. 
This corresponds to approximately 10% of particles having scat- 
tered over a Hubble time at this radius. However, the changes are 
small compared to the change where V x (10 Gyr) > 1. We note 
here that the most massive halos in Dave et al. (2001 ) also seem 
to show the same qualitative behavior. However, even for large val- 
ues of the deviation from sphericity is significant, a fact that 
is in some disagreement with the simulation results of Dave et al. 
( |2001 ) > where c/a > 0.9 for their most massive halo. We specu- 
late that part of this could be due to the differences in the way the 
ellipticity was estimated and part could be due to the smaller box 
run by Dave et al. (2001 1, which implies a quieter merger history. 
It is well known that CDM halos have anisotropic velocity ellip- 
soids and elongated shapes that are driven partially by directional 
merging (Allgood et al. 2006). These mergers provide a source of 
anisotropy that needs to be overcome by scattering in order for ha- 
los to reach sphericity. We also note that the energy transfer facili- 
tated by self-interactions would lead to an isotropic velocity disper- 
sion tensor and that does not necessarily imply a rounder halo. To 
make this connection between isotropic velocity dispersion tensor 
and a rounder halo, previous analytic estimates have relied on the 
simulations of Dave et al.| ( [200T| >. 

The difference may also be a numerical artifact: we find that 
halo shapes only converge if there are at least 10 4 particles in the 



halo and only for radii r > 2r r e i ax (see Se c.[2]l. Most of the ha 



Dave et al. 



(2001) only have W 5 



los used for shape estimates in 
particles in the virial radius. Their largest halo does have more 
than enough particles for reliable shape estimates. However, for 
this massive halo, |Dave et al.| j2001) show shape measurements 
at radii much smaller compared to the convergence radius than we 
do. In our simulations, we find that halos appear artificially round 
below the convergence radius. Our shape measurements are con- 
sistent with Dave et al. (2001)'s massive halo for radii above the 
convergence radius. 

The other major effect of energy transfers due to scattering 
is to create a core (see Rocha et al. (2012)) and deep inside a con- 
stant density core, we must have c/a — > 1. Hence we expect the log 
slope of the density profile to correlate strongly with the halo shape. 
Thus, instead of the density profile scaling as p ~ r _1 in the inte- 



rior as expected for CDM (Navarro, Frenk & White 1997; Navarro 
|et al.|2004[|2010} , the density profile of SIDM halos plateaus close 
to the center. We use the following proxy for the negative log slope 
of the density profile: 



7 = 3- 47rp(r)r 3 /M(r), 



(3) 



where M(r) is the mass enclosed within radius r and 7 = 
— d log p/d log r for a power-law density profile. Based on the pre- 
vious figures, we expect halo shapes to become increasingly round 
as 7 < 1, and that the CDM and a/m — 0.1 cm 2 /g halos should 
not dip below 7 « 1. In Fig. [6] we show c/a as a function of the log 
slope of the density profile, as approximated using Eq. {3}, for host 
halos in our highest mass bin. We obtain similar results for halos 
of smaller mass, but only show the highest mass bin because these 
halos are the best resolved. 

We find that, indeed, SIDMi halos become significantly 
rounder as 7 < 1 and become almost completely round when 7 
gets much smaller than 0.5. Interestingly, we also see that c/a de- 
viates strongly from CDM even for relatively large values for the 
log slope, in regions of the halo in which the scattering is not effi- 
cient at changing the radial density profile. This is a consequence 
of the fact, as shown in Fig. [5] that it does not take a lot of scatters 
to start rounding out the halos, although it takes multiple scatters 
for the halos to acquire c/a axis ratios in excess 0.8. Thus, the ef- 
fects of scattering are apparent for a/m — 1 cm 2 /g even when 
p ~ r~ 2 and the density profile is unaffected by scatterings. How- 
ever, the observational importance of this behavior is mitigated by 
two factors - the change for 7 > 1 is mild and easily within the 
scatter in ellipticities seen in CDM. 

In summary, we find that it only takes a modest local scat- 
tering rate per particle, T(r) > O.lHo, to start changing the 
three-dimensional halo shape within radius r with respect to CDM. 
We find that SIDMi halos gets significantly rounder compared to 
CDM predictions when the negative log-slope of the density pro- 
file 7 <C 1 in regions where dark matter particles (on average) have 
had at least interactions in a Hubble time (Tx (10 Gyr) > 1) . Our 
results show that even in the limit of one or more scatterings, the 
halo shapes retain some of their initial triaxiality. 



4 COMPARISONS TO OBSERVATIONS 
4.1 Defining Observables 

There are a number of ways of quantifying deviations of mass dis- 
tributions from spherical or axial symmetry. In the previous sec- 
tion, we quantified the deviations in terms of c/a, the ratio of 
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Figure 5. Median axis ratio c/a for host halos as a function of the local scattering rate modulo the cross section: pv rrnB ~ Y{a/m)~ . Smaller values 
correspond to the outer halo, where the density and scattering rate are low. The quantity is scaled by 10 Gyr cm 2 /g, such that 1 in these units means that each 
particle has roughly one interaction per 10 Gyr in SIDMi (blue dashed line) and one interaction per 100 Gyr in SIDMo.i (green dotted line). The black solid 
line is CDM. 



the semi-minor to semi-major axes determined from the modified 
three-dimensional moment-of-inertia tensor, Eq. This is rarely 
a practical shape estimate observationally. Instead, there is a more 
suitable measure of halo ellipticity or triaxiality for each type of ob- 
servation. The relationship between these measures is non-trivial, 
so care must be taken to compare theory to observation appropri- 
ately. From paper to paper, the definition of "ellipticity" can change 
significantly. In order to facilitate these comparisons, we define 
three distinct measures of asymmetry in this section and go on to 
use them for specific comparisons to observational studies in Sec. 
|4.2| |4.3| and |4.4| The symbols we use for these shape definitions 
are summarized in Table 2. 

For strong lensing, which we discuss in Secs. |4.2| and |4.4| what 
matters is the deviation of the convergence from axial symmetry. 
The convergence is k — E/E cr , where E cr is the critical surface 
density for creating multiple images. For an axially symmetric sys- 
tem the convergence will be k(6), and depend only on the angular 
distance on the sky 9 from then lens center. Once axial symmetry is 
broken, one must consider k(6, (f>), where cj> is the azimuthal angle 
that rotates on the sky. Miralda-Escude ( 2002 1 used the following 
quadrupole approximation to fit the surface density of MS 2137-23, 



K(6,<t>) = K (9) 



(4) 



where Ko(0) is the convergence averaged along azimuthal angle 
and e quantifies the amplitude of deviation of the convergence from 
axial symmetry. Since the normalization of k and the angle 9 de- 
pend on the source, lens, and source-to-lens distances, which gener- 
ically vary, we quantify deviations from axial symmetry in terms of 
surface density E(ii, <f>), where R is the two-dimensional physi- 
cal radius in projection. Using Eq. |4}, we define the measure of 
ellipticity 



e'(R, cj>) = 



2(E o (fi)-E(fi,0)) 
dEo/di? 



(5) 



which should be equivalent to e cos(2<j!>) if the quadrupole expan- 
sion of the two-dimensional surface density is approximately cor- 
rect. Here, En(ii) is the azimuthally averaged surface density. 

The second type of measure used to quantify deviations from 



spherical or axial symmetry in lensing arises from the extension of 
the double pseudo-isothermal sphere, 



p(r) 



Po 



(6) 



to allow for deviations of the surface-mass density from axial sym- 
metry (double pseudo-isothermal elliptical, or dPIE: see Richard 
|et al. |2010) , 



2G r, 



fcut 
cut f'core 



with 



f 2 = [(2-e)( a 



,)/2] 2 + [(2- e )( J/ -y c )/(2-2e)] 2 



(8) 



The dPIE profile has the properties that p ~ const for r <C r core , 
p oc r~ 2 for r C ore < r -C r c uu and p oc r~ 4 for r » r cut , 
so that the total mass is finite. The center of the halo in projection 
is denoted by (x c , y c ), with the 2>direction aligned with the major 
axis of the distribution. This surface-density profile is often used in 
fits of the shapes of galaxies or dark-matter halos in clusters that 
strongly lens background galaxies. 

The third way we will quantify halo shapes observationally is 
in terms of similar spheroids (see, e.g., Section 2.5 of |Binney &] 
Tremaine 2008 ). For spheroids, we can define an ellipticity param- 
eter 



e = 1 



(9) 



where b is the semi-minor and a the semi-major axis. If the z-axis 
is the symmetry axis, the semi-major axis is given by 



2 D 2 . Z 

a — H + 



(l-c)*' 
for oblate spheroids, and 

R 2 



oblate 



+ z , prolate 



(10) 



(11) 



(1 - e )a 

for prolate spheroids. Similar spheroids are those for which the 
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Table 2: Summary of shape definitions used in observational comparisons 



Symbol 


Defining Equation 


Description 


Relevant Figure 


e'(R, 0) 




Deviation from axial symmetry in leasing convergence maps 


Fig.g 


e 


Eq|8] 


Ellipticity in dPIE surface density fits to lensing signal 


Fig.[ll] 


c 


Eqf] 


Ellipticity of similar spheroid used in X-ray studies 


Kg.§ 
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Figure 6. Median halo shapes as a function of the estimated log slope of 
the radial density profile. Lines have the same meaning as in Fig. [3] 



density profile may be described in terms of p(a) and e is fixed 
throughout the body. The isopotential surfaces of such spheroids 
are rounder at large distances if the body is more centrally concen- 
trated (see the discussion in Sec. 2.5 of Binney & Tremain e[2008| ). 

As summarized in Table 2, e' is the surface-density shape def- 
inition we use in Sec. |4.2| e is the X-ray motivated shape definition 
we use in Sec. |4.3| and e is the lensing-fit shape definition used in 
Sec.l44l 



4.2 Revisiting Mira lda-EscudilpOO^ 

Galaxy clusters are great places to look for the effects of velocity- 
independent SIDM because one may typically achieve much higher 
values of p(r)v rma (r) for fixed r/r v i T . In addition, there are many 
different probes of the mass distribution of clusters that span an 
enormous dynamic range of radial scale — stellar kinematics and 
strong lensing towards the center of the cluster, weak lensing and 
X-ray gas distributions throughout the halo volume, and weak + 
strong lensing maps of the matter distribution around individual 
galaxies in the cluster ( San d et al.|2008||Newman et al.|2009||201 1| 



|Kneib & Natarajan|201lT 7¥ is no surprise the tightest constraints 
on velocity-independent SIDM emerged from cluster studies, and 
a revisit of the tightest of these constraints is the subject of this 
section. 

The strongest published constraint on velocity-independe nt 
SIDM, a/m < 0.02 cm 2 /g came from 



Miralda-Escude 



ipendent 
j2002| 's 



mated virial mass of ~ 8 x 10 14 M Q {Gavazzi|2005| , and its mass 
distribution has also been studied by|Fort et al. 1 1992|>,|Mellier, Fort] 
& Kneib|(T993) >, |Miralda-Escude| ( fT995^|Gavazzi et al.| i |2003| >, and 
Sand et al.|(;2008 1. There are two strongly lensed galaxies that pro- 



study of the galaxy cluster MS 2137-23. This cluster has an esti- 



duce a total of five distinct images: one source has a radial image 
at 9 ~ 5" from the center of the brightest cluster galaxy and an 
arclet at 8 = 22.5". The other source has a large tangential arc and 
two arclets all at about 6 — 15" from the brightest cluster galaxy 
center, which corresponds to 70 kpc. In order to reproduce both the 
relative magnifications and the alignments of the images in the sky, 
the surface density must deviate from axial symmetry at 70 kpc. 
Quantitatively, it means that the parameter e in Eq. which cor- 
responds to the amplitude of e' given in Eq. {5), must be e ~ 0.2 
at R = 70 kpc. This figure is largely driven by the tangential arcs 
and associated arclets (Miralda-Escude 1995) 

Based on the ellipticity at 70 kpc and an argument that the 
dark-matter surface density should be approximately axial for a 
typical particle collision rate T > Hp , |Miralda-E scude (2002} as- 
serts that F(70 kpc) < Hq. Using the fact that the tangential arc 
should lie at approximately where the mean interior convergence 
R — 1 (or an estimated critical density E cr = 1 g/cm 2 ) and the 
rough approximation p(r) ~ S(r) /r, [Miralda-Escude| estimates 
the three-dimensional density p(70kpc). Using the velocity dis- 
persion of the brightest cluster galaxy at the center of the halo as 
a proxy for u rms , Miralda-Escude uses Eq. |TJ to determine a limit 
on a/m, which is found to be cr/m < 0.02 cm 2 /g. 

We get a sense that this line of reasoning may be flawed when 
we examine the surface-density plots of one of our most massive 
halos in Figs. [T] and [2] and our findings of Sec. [5] First, recall that 
our results show (cf., Fig. [5} that the inner halo shape retains some 
triaxiality even when F > Ho- Second, the surface density includes 
all matter along the line of sight, not just the material within r < R. 
Thus, the surface density at small R/r v i T includes a lot of material 
with large r, far out in the halo where SIDM scatters are unimpor- 
tant. This material is still quite triaxial. Moreover, SIDM also cre- 
ates cores, which means that the outskirts of the halo have an even 
greater weight in the total surface density than if the halo were still 
cuspy at the center. Empirically, we see that the simulated surface 
densities in Figs. [TJ and [2] are quite elliptical. This point becomes 
more and more important as the size of the core becomes smaller. 
All of these things suggest that the constraint reported in Miralda- 
Escude ( 2002 1 is far too high. 

When attempting to quantify the constraints on SIDM from 
MS 2137-23 using our simulations, we run into the following 
problem: the largest halo in our simulations has a virial mass 
M v i r = 2.2 x 1O 14 M0, a factor of approximately four smaller 
than the estimated virial mass of MS 2137-23. Moreover, we do 
not know the orientation of the principal axes of the cluster with 
respect to the line of sight. In order to make the comparison, we 
do two things. First, we can use virial scaling relations to estimate 
the radius at which p(r)v rma (r), a proxy for the scattering rate (see 
Fig. |5j, has the same value as it would for a radius of 70 kpc in 
MS 2137-23, in other words, we look for the radius at which the 
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Figure 7. The surface-density ellipticity e'(R, <f>) as denned in Eq.j5]for the same halo as in Figs. [T] and [2] Shown are projections along the three principle 
axes, evaluated at R = 35 kpc and plotted as a function of azimuthal angle. As in the figures in Sec. [3] the solid black line corresponds to CDM, the blue 
dashed line to cr/m = 1 cm 2 /g, and the green dotted line to a/m = 0.1 cm 2 /g. The dot-dashed red line shows the minimum amplitude of e'(4>) required in 
order to explain the image positions of MS 2137-23. We do not show the lines for the a/m = 0.03 cm 2 /g simulation because it is indistinguishable from the 
CDM line. The angle (f> is defined such that = corresponds to the longer axis of the surface-density distribution. 



SIDM scattering rate should be comparable to the radius at which 
|Miralda-Escude1 $2W2\ finds the constraint for MS 2137-23. For 
our M vir ~ (1 - 2) x 1O 14 M halos, r = 35 kpc is roughly the 
point at which the scattering rate is similar to that at 70 kpc in MS 
2137-23. This is outside r m i n for these halos, so we trust the shape 
measurements. Second, we look at several projections of the ha- 
los. We calculate E (-R), <j>), and hence e'(cj>) for the various 
projections of the halos. 

We show an example of e'{4>) curves for lines-of-sight along 
the principal axes of the halo moment-of-inertia tensors of one of 
our largest halos in Fig [7] the same halo shown in Figs. [T] and [2] 
As in previous figures, the solid black line denotes the CDM re- 
sult, the blue dashed line the result for a/m = 1 cm 2 /g, and the 
green dotted line the result for a/m — 0.1 cm 2 /g. We do not show 
the curve for a/m — 0.03 cm 2 /g even though it is the cross sec- 
tion closest to the Miralda-Escude (2002) constraint because it is 
indistinguishable from the CDM line. The red dotted line shows 
the minimum amplitude of e'(<f>) required for the lens model of MS 
2137-23. The e'(e6) curves of the other massive halos look simi- 
lar to the curves for the halo shown in Fig. [7] What we find is that 
even the a/m = 1 cm 2 /g curve generally satisfies the MS 2137-23 
constraint. Therefore, we find that the Miralda-Escude (2002) con- 
straint is in fact overly constraining by two orders of magnitude. 
While we do not simulate cosmologies with cr/m > 1 cm 2 /g, and 
thus cannot set a quantitative upper limit on the SIDM cross sec- 
tion, we may conclude that a/m = 1 cm 2 / g is not ruled out by 
MS 2137-23. 

There are a few caveats to this conclusion, which we do not 
believe will significantly alter our claim. First, none of our sim- 
ulated clusters are as massive as MS 2137-23. This precludes us 
from doing a detailed comparison of the projected densities in 
a/m = 1 cm 2 /g to that required to explain the arcs in MS 2137- 
23. It is, however, informative to do a simple calculation to gauge 
the importance of this effect. We appeal to the results in the com- 
panion paper (Rocha et. al. 2012) that show that the density profile 
in SIDM for a/m — 1 cm 2 /g is well fit by a Burkert profile ( Burk- 



|ert|1995) , and that this profile deviates from CDM density profiles 
at radii smaller than half the Burkert scale radius. At this point, 
rb/2 or 0.35(r ma x/21.6 kpc)~ 0,08 , the density profile becomes al- 
most constant. We use the Vmax — fmax relation seen in our CDM 
simulations and the NFW profile to compute the projected mass 
within 70 kpc as the CDM prediction and compare that to the SIDM 
prediction by computing the same projected density but assuming 
fhatVr < rb/2, psiDM(r) = pcDM(n>/2). This computation re- 
veals that the projected density in a/m = 1 cm 2 /g model should 
be about 30% lower than in CDM for the median (1 - 5) x 10 15 M Q 
halos. Even for the CDM case, however, the projected density is 
about a factor of 2 smaller than the estimated E cr = 1 g/cm 2 . 
Clearly, the estimated projected density (as opposed to the shape) 
could be a significant constraint on the a/m — 1 cm 2 /g model. 
This simple computation motivates further work along these lines 
with more realistic SIDM density profiles, inclusion of scatter and 
the uncertainty in the halo virial mass. 

We have started much larger-scale simulations in order to 
study clusters in more detail, which will include an investigation 
of strong-lensing cross sections as well as ellipticity. In terms of 
the ellipticity function e' (</>), it is not clear in which direction our 
results will go for simulated 8 x 1O 14 M0 halos. On one hand, large 
halos are more triaxial than small halos, so we would expect that if 
anything, we would be underestimating the degree of ellipticity in 
the surface-density distribution. On the other hand, large halos have 
larger and rounder cores for velocity-independent SIDM compared 
to their lower-mass cousins. This might drive the constraints the 
other way, although a larger core implies that the outskirts of the 
halo get weighted more heavily along the line-of-sight integral of 
the density for the surface-density calculation. Finally, several au- 
thors have noted that MS 2137-23 is actually unusually round for a 
galaxy cluster (Gavazzi 2005, Sand et al. 2008). And intriguingly, 
the modeling of both |Miralda-Escude| ( |1995fr and |Sand et af| ( 2008 ) 
indicate that a cored dark-matter radial density profile is preferred 
over a strongly, CDM-like cuspy profile for the dark-matter halo. 
Such a cored profile is more in line with what we would expect 
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Figure 8. Distribution of halo ellipticity e as defined in Eq.[9]from fits 
for halos with M = (3 - 10) X 10 12 Af Q . The black histograms indi- 
cate ellipticities for CDM halos, and the blue shaded histograms show the 
ellipticities for halos with a/m = 1 cm 2 /g, and green dashed histogram 



shows ellipticities if a/m = 0.1 cm /g. 

shows the center and width of the best-fit e and uncertainty from 



The dark blue line with arrow 
Buote| 



|et al. 1 (2002 ). The shapes are found in using the Dubins ki & Carlberg| < 



19911 



weighted moment-of-inertia method described in Sec.[2]and approximated 
as oblate or prolate spheroids. The initial minimum radius of the shell for 
the shape measurement is set to r m ; n = 8.5 kpc, and the maximum radius 
to r ma x = 14 kpc. See text for details. 



for SIDM. Tighter constraints on a/m, or a measurement should it 
be non-zero, will come from an ensemble of clusters, including the 
most triaxial of them, a point to which we return in Sec.|4.4| 



4.3 Shapes from X-ray observations of elliptical galaxies 



Constraints on dark-matter scattering have also been made us- 
ing the shapes of significantly smaller dark-matter halos, M v j r ^ 
M ra jFeng et al.|2009[|Feng, kaplinghat & Yu|2010 



10 



10 1 



|Buckley & Fox|2010| |McDermott, Yu & Zurek 201 1 1. In this case, 
the observations consist of X-ray studies of the hot gas halos of 
elliptical galaxies (Buote et al.|2002| l. 

One may use the shape and twisting of the X-ray isophotes 
to learn about the three-dimensional shape of the dark-matter halo. 
Interestingly, halo shapes may be constrained with imaging data 
alone and do not rely on temperature-profile modeling, which is 
required for mass-density determinations. This "geometric argu- 
ment" was first made by Binney & Strimpel (1978) and subse- 
quently applied by a number of authors in the study of ellipti- 
cal galaxies and galaxy clusters (Fabricant, Rybicki & Gorenstein 
[T984l[Buote & Canizares|1994||1996||1998b||Buote et al.|2002| >. 
The geometric argument is the following: for a single-phase gas in 
hydrostatic equilibrium, the three-dimensional surfaces in the gas 
temperature T, gas pressure p g , gas density p g and gravitational po- 
tential $ all have the same shape (Binney & Strimpel|1978| Buote| 



& Canizares 



1998a 1. Surfaces of constant emissivity jx oc p g have 



the same three-dimensional shape as the isopotential surface. 

On the other hand, if the spectral data are available (as they 



typically are with the Chandra telescope), then the temperature- 
profile data may be used with the assumption of hydrostatic equi- 
librium to fit the X-ray isophotes. |Buote e~ al. (20021 use both ap- 
proaches to model NGC 720. 

X-ray isophotes are a good probe of the shape of the lo- 
cal matter distribution for the following reasons. First, the total 
mass profile (galaxy + gas + dark- matter halo) is nearly isothermal 
(p 



n 1 1 1 r 

r ) for elliptical galaxies (Humphrey et al. 2006; Gavazzi 



|et al.|2 007 ). This means that the shape of the isopotential contours 
reflects the shape of the matter distribution at the same radius, al- 
though typically the isopotential surfaces are significantly rounder 
than the isodensity contours (see, e.g., Binney & Tremaine 2008| >. 
Second, the density profile of gas tends to be pretty sharply cuspy 



'). Since the emissivity goes as jx oc p g , and the mor- 



phology of jx traces that of the gravitational potential (see the dis- 
cussion on the geometric argument above), this means that most of 
the X-ray emission along the line of sight is concentrated at radii 
close to the projected radius. Thus, the X-ray isophotes indicate 
the shape of the matter distribution at radii similar to the projected 
radius. 

In this section, we focus on the shape measurement of the 
X-ray emission about the elliptical galaxy NGC 720. The inferred 
shape of the matter distribution of this system has been used to set 
constraints on self-interacting dark matter in the recent past (|Feng| 



et~aLp 009; Feng, Kapling hat & Yu|2010| [Buckley & Fox||2010[ 



McDe rmott, Yu & Zurek|201 1| >. The data set used by |Buote et al.] 
1 2002) for the shape measurement of NGC 720 is a 40 ks expo- 
sure of the inner 5' of the galaxy with the ACIS-S3 camera on 
the Chandra telescope. The data included in the fit are contained 
within a 35" — 185" annulus from the center of the galaxy, which 
corresponds to « 4.5 — 22.4 kpc. |Buote et al. (2002) estimate the 
three-dimensional isopotential shapes in the following way. First, 
they investigate a mass-follows-optical light (M oc L«) mass pro- 
file for the gravitating mass using the geometric argument. They 
use a spheroidal Hernquist ( 19901 model for the stellar mass of the 
galaxy, with structural parameters determined by deprojecting the 
optical image to three dimensions. This is to test if the stars in the 
galaxy may be sufficient to provide a gravitational potential for the 
gas. They find that the isophotes are rounder than observed beyond 
about R e (ss 50" or 5 kpc). This suggests that there must be sig- 
nificant ellipsoidally distributed mass extending well beyond the 
effective radius of stars, since isopotential contours become round 
as the distance to the main ellipsoidal mass profile increases. 

Next, |Buote et ak] (2002) use the fact that they find the tem- 
perature profile of the halo gas to be isothermal to find the three- 
dimensional X-ray emissivity distribution (jx oc pi) directly from 
the equation of hydrostatic equilibrium: 



Pg(r) 



— exp 



/xm p <J>(0) 
k B T 



(#(r)/«(0) - 1) 



(12) 



where /1 is the molecular weight of the gas atoms and m p is the 
proton mass. They used three similar spheroid models for the mass 
distribution of the galaxy (baryons plus dark matter) to find the 
gravitational potential 4>(r). The three models were pseud oisothe r- 
mal (p(a) oc 1/(1 + (a/a p ) 2 ) where a is defined in Eqs. jlQ|llf ), 
Navarro-Frenk- White (Navarro," Frenk & Whitel[T997| and Hern- 
quist (Hernquist 1990). They fit both oblate and prolate spheroids, 
and assumed that the symmetry axis was in the plane of the sky, as 
this leads to the most elliptical isophotes for a given spheroid. Once 
the potential was found, they calculated the emissivity distribution 
in three dimensions and integrated along the line of sight. A \ 2 
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statistic was used to test the quality of the model fits to the X-ray 
surface-brightness map. 

The best-fit models were the pseudoisothermal ones, with the 
oblate and prolate spheroids having nearly identical goodness-of- 
fit. The best-fit ellipticity e = 0.37 ± 0.03 for the oblate spheroid, 
and e = 0.36 ± 0.02 for the prolate spheroid. The flatness of the 
ellipticity as a function of radius in the region of interest drives the 
best-fit density-profile for isofhermality. Only for a similar isother- 
mal spheroid are the X-ray isophotes so constant in ellipticity as- 
suming density profiles that are similar spheroids. 

We compare our simulations to the B uote et aL] ([2002 ) results 
by using the weighted moment-of-inertia tensor ellipsoidal shape 
estimator of Sec. [2] This method weights particles in the region 
of interest equally in estimating the shape of the mass distribu- 
tion, and does not preferentially weight particles near the edge as 
the true moment-of-inertia tensor does. The initial region in which 
the ellipsoidal shapes are estimated is a spherical shell of radius 
r-min < r < r- max , where r m in is our numerical limit on the shape 
convergence radius and r max = 14 kpc. We choose this value of 
fmax as a compromise between finding the ellipticity at small radii, 
where differences with CDM are highest, and having enough par- 
ticles in the region of interest for a robust and unbiased shape es- 
timate. The shape of the region of interest is deformed in each it- 
eration of the weighted moment-of-inertia calculation to reflect the 
major axes of the morphology of the tensor at that iteration, keeping 
the morphology of the region of interest ellipsoidal and fixing its 
volume. This region is within the core radii of the a/m = 1 cm / g 
halos, and is our best approximation to the shape of the innermost 
part of the dark-matter halos, the parts relevant to the study of NGC 
720. For reference, the 10 12 , 1O 13 M0 halos have median core radii 
16, 43 kpc, respectively (Rocha et. aL, 2012). 

In order to find the spheroidal ellipticity (described in Eq. (|9]l) 
from these ellipsoidal fits, we use the ellipsoidal axis ratios to de- 
cide if the halo is oblate or prolate. If prolate, we take axis ratio 



of the spheroid 1 — e, cf. Eqs. 



IOi & 
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to be y/bc/a, i.e., the 



spheroidal semi-minor axis is the geometric mean of the ellipsoidal 
minor and intermediate axes. If oblate, we set 1 — e = c/yab, i.e., 
the semi-major axis of the oblate spheroid is set to be the geomet- 
ric mean of the major and intermediate axes of the ellipsoid. Both 
these choices set the spheroidal volume equal to the ellipsoidal vol- 
ume. A more realistic analysis would compute the integrated jx 
directly from the potential of the simulated halo and compare that 
to the observations. However, there are significant other uncertain- 
ties (discussed below) that argue against such an approach being 
more fruitful. 

In Fig. [8] we show the ellipticity distribution for all halos in 
our 25 /i _1 Mpc boxes for CDM and SIDMi that have a virial mass 
within the l a uncertainties of the mass modeling in Humphrey 
i2006] for NGC 720, Af vir 



et al. 



(6.6:! 



x 10 12 M Q . Note 



that we do not weigh the e distribution by the error distribution for 
the virial mass presented in Humphrey et al. ( 2006 1. However, the e 
distribution in our simulated halos is not a strong function of virial 
mass in this mass range. We show the CDM e distributions with 
the black histograms, the a/m = 1 cm 2 /g e distributions with the 
cyan shaded histograms, and a/m = 0.1 cm 2 /g with dashed green 
histograms. We find that the inferred ellipticities are approximately 
independent of r max out to approximately 25 kpc, with there being 
a small tail at higher e for the a/m — 1 cm 2 /g halos since the 
core radii are in the 20-40 kpc range and halo shapes are less af- 
fected by scatterings at radii larger than roughly the core radius. For 
a/m — 0.1 cm 2 /g, the core sizes are of order r m i n or smaller, so 

the shape measurements are relatively insensitive to r max range. A 1 2006 >. The core sizes for a/m = 0.1 cm 2 /g are smaller, ~ 7kpc 



quick review of Fig.[3]also shows that the relative independence of 
the ellipticity estimate out to 0.1r v i r is to be expected in each cos- 
mology. For the a/m — 1 cm 2 /g halos, we exclude those halos for 
which either poor centering of the halo or ongoing merging makes 
the halos appear artificially flattened. For the other cosmologies, 
the relative cuspiness of the central density leads to more accurate 
halo centering. 

Although there are significant differences in the ellipticity 
distributions of the CDM and SIDM halos, the observed elliptic- 
ity (downward arrow) is just within the ellipticity distribution for 
a/m — 1 cm 2 /g. The fact that this magnitude of ellipticity is still 
in the distribution for a/m — 1 cm 2 /g is also apparent from the 
middle panel of Fig. [3] which shows that the three-dimensional axis 
ratios at r/r v - a ~ 0.04 (approximately the radius corresponding to 
the 2D region of interest for NGC 720) can accommodate signif- 
icant ellipticity. Note that radius is also close to the typical core 
radius where there we expect roughly 1 interaction per particle on 
average. This implies that NGC 720 could be an extreme outlier in 
SIDMi model and thereby consistent with observations. 

However, there are some issues with this interpretation. One 
issue is that the ellipticity measured for the dark matter halo seems 
to be constant down to roughly ~ 5 kpc in the data, which is 
about 1% of the virial radius for the mean virial mass for NGC 
720 ( [Humphrey et al.[ 2006 1. The middle panel of Fig. [3] shows a 
sharp rise in axis ratio when the number of interactions gets above 
unity and hence our results are likely to change if we were able push 
deeper into the core with higher-resolution simulations. A second 
issue is that the measured halo mass for dark matter within 10 kpc 
((Humphrey et al. 2006 1 is a factor of 2-3 larger than even the me- 
dian mass in our CDM simulations for halos with virial mass in the 
range preferred by Humphrey et al. (2006) fits. Thus we should also 
cut the histogram in Fig.[8]based on, say, the mass within 10 kpc (or 
some other measure of the concentration). These halos, by virtue of 
their larger densities, will also be rounder. The discrepancy in av- 
erage dark-matter density within 10 kpc may also be a sign that the 
inner part of the dark-matter halo has been compressed as a conse- 
quence of assembly of the elliptical galaxy, a compression that has 
been observed in other studies of elliptical galaxies ( jSchulz, Man-| 
|delbaum &'P admanab han|2010[|Dutton et al.|201 1| >. Moreover the 
presence of the stars and gas increase the velocity dispersion of the 
dark matter in the central parts and increase the scattering rate - an 
effect that is not captured in our SIDM simulations. 

The resolution of these issues lies in a better comparison to the 
data, which in turn requires a bigger box higher-resolution simula- 
tion to probe deeper into the halo and gain more statistics. It will be 
important to include baryons to see how their presence may affect 
halo properties. With such a halo catalog in hand, it will be inter- 
esting to do a more careful comparison to the X-ray data of NGC 
720 and other large nearby ellipticals. The addition of other ellip- 
ticals would be crucial. With only one object, the spread we see in 
ellipticities may be hard to overcome (although it may be smaller 
if we also cut on concentration as discussed above). If we had an 
ensemble of shape measurements, we would be able to set tighter 
limits on the SIDM cross section. 

While our comparison to |Buote et al.| ( |200"2| l is not sufficiently 
sharp, the weight of the arguments suggests that a/m = 1 cm 2 /g 
is not likely to be consistent with the measured shape of NGC 720. 
However, based on our existing simulations, a/m = 0.1 cm 2 /g is 
as consistent as CDM for the shape of the NGC 720 isophotes. It is 
interesting to note in this regard that there is no hint of a large core 
(~ 30 kpc) in the results of Buote et al. ( 2002) or Humphre y et al.| 
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for the same virial mass range (Rocha et al, 2012), comparable to 
the effective radius of the stars in NGC 720. Thus the inferred dark 
matter density profile in NGC 720 may be a better way to search 
for effects of self-interactions. 

To amplify the point about the dark matter density further, we 
note that the median central (maximum) density for SIDM with 



Humphrey et al.H2006) 



a/m = 0.1 cm 2 /gis O^Mo/pc 3 , while|] 

infer an average density of O.O4M0/pc 3 within 10 kpc. At 5 kpc, 
the inferred average density is 0.1Mo/pc 3 , still within a factor of 
two (expected from scatter) of the predictions. This lends credence 
to the argument that an analysis focused on SIDM predictions of 
an ensemble of nearby X-ray-detected elliptical galaxies could be 
a fruitful way to look for signatures of or constrain SIDM. It is 
also worth noting that the shape distribution of a/m — 0.1 cm 2 /g 
is visibly different from CDM and perhaps an ensemble of X-ray- 
shape measurements could resolve the differences. 

Our analysis argues for the conclusion that SIDM cross sec- 
tions with a/m > 0.1 cm 2 /g can be hidden in X-ray data, and 
that previous constraints on SIDM using these data are overly strin- 
gent. Feng, Kaplinghat & Yu (2010 ) assumed that the Buo te et alT] 
( |2002| ellipticities described the halo shape at the inner radius of 
R = 4.5 kpc and that F ~ Hq was required to make the halo 
spherical as indicated by the results of Da v6~et al.H206"T) . Our re- 
sults indicate that this interpretation is flawed because (a) SIDM 
halos retain significant triaxiality in the region where T ~ Ho and 
(b) the scatter in SIDM halo ellipticities is large. In order to use an- 
alytic arguments to constrain self-interacting dark matter models, 
they should be tuned to reproduce the the distribution of axis ratios 
seen in simulations. 



4.4 The future of cluster lensing constraints 

In the future, far better constraints on SIDM will come from statis- 
tical studies of galaxy-cluster lens samples rather than the analysis 
of individual objects. In this section, we focus on statistical studies 
of the shapes of relaxed clusters. 

There are a number of ongoing and future observational pro- 
grams designed to characterize the mass function of and mass dis- 
tribution within galaxy clusters (e.g., 



tions 



2010 



LSST Science Collabora- 
[2009) |Gill et al.||2009"l |Plagge et al.||2010| |Richard et al 



Planc k^ollaboration et al.|201 la|b|c| Pillepich, Porciani & 



Reiprich||2012| IMarriage et al.||2011||Viana et al.|2012| |Postman 



et al.||2012) . Modulo the effects of baryons, a smoking-gun sign 
of SIDM would be for the mass function of galaxy clusters to look 
identical to CDM, but with lower mass density and rounder surface- 
mass distributions at the centers of the clusters. One would use 
the ensemble of galaxy-cluster data to compare with simulations 
of clusters in CDM and SIDM (with various elastic scattering cross 
sections). 

In this study, we will consider shape-based constraints from 
the initial results of the Local Cluster Substructure Survey (Lo- 
CuSS; PI: G. Smith), a multi-wavelength follow-up program of 
165 low-redshift clusters selected from the ROSAT All-sky Sur- 
vey catalog (Ebeling et al. 2000; Bohringer et al. 2004; Richard 
|et al.|2010|. Twenty cluster s were used for the first mass-modeling 
study {Richard et al.|2010[ >. These were selected because they had 
been observed with the Hubble Space Telescope (HST), could be 
followed up spectroscopically at Keck, and were confirmed to have 
strongly lensed background galaxies. The details of the mass mod- 
eling of these clusters are presented in Sec. 3 of Richard et al. 
(2010). For our purposes, the key fact is that the surface density of 
the cluster-mass dark-matter-halo component of the lens model was 



parametrized in terms of the dPIE profile, Eq. |7|. This study should 
be taken as an example of the power of using statistical studies of 
clusters for SIDM constraints, with an emphasis on the constraints 
on halo shapes. 

In order to compare the LoCuSS clusters to simulations, we fit 
the surface densities of our most massive clusters in the 50 h ~ 1 Mpc 
boxes for our CDM, SIDMo.i, and SIDMi runs using the dPIE 
surface-density profile, Eq. |7|. We fit the surface densities of the 
halos projected along the principal axes of the moment-of-inertia 
tensors of the mass within the virial radius, and perform the fits 
within annuli R m in < R < i? m ax- The inner radius of the an- 
nuli is chosen to be _R m i n = 20kpc//i, since this is the three- 
dimensional shape convergence radius (see Sec.[2](. We set the outer 
radius i? ma x = 50 kpc/ h since most of the lensing arcs in Lo- 
CuSS are in the range 20 kpc < R < 100 kpc if projected into 
the plane of the sky at the position of the galaxy clusters. We fix 
r C ut = 1000 kpc, which is what is typically done with the LoCuSS 
clusters. In practice, the parameter constraints are insensitive to this 
choice, as all the data are well within this radius if projected onto 
the sky. The free parameters that we fit are: ao, the characteristic 
velocity dispersion of the cluster; r core , the pseudoisothermal core 
radius; e, the cluster ellipticity (as defined in Eq. {8j); and the posi- 
tion angle 6. We fix the center of the cluster to that inferred by AHF 
The surface-density parameters were fit using the downhill simplex 
algorithm in the scipy python module using the likelihood 



C= P(N obs \N{a ,r c 



,o)) n a 



(13) 



where P(N bs\N(ao,r core ,e,9)) is the Poisson probability of 
finding N t s simulation particles within the annulus given a dPIE 
model with parameters (ao, r cor(! , e, 8), for which TV simulation 
particles are expected, d is the probability of finding a simulation 
particle i at position (x, y) given those same parameters, 



C 



E(x,y\a ,r core ,e,6) 



(14) 



f , dxdyY,(x' y') ' 

J annulus ° * ' •> 1 

where T,(x, y) is defined in Eq. |7|. 

Examples of the dPIE fits are shown in Figs. [9] and [TO] in 
which we show the surface densities of one massive halo (A/ V ir = 
1.8 x 1O 14 M0) along the major and intermediate principal axes 
of the halo moment-of-inertia tensor as well as the best-fit dPIE 
surface densities. The central regions of the halos are masked for 
R < fmin, the region where the shape profiles in three dimensions 
are not converged. The halo appears rounder and more dense if pro- 
jected along the major axis rather than the intermediate axis, just as 
we saw for another halo in Figs.[T]and|2] While there are noticeable 
differences between the CDM and a/m — 1 cm 2 / g surface den- 
sities at small projected radii, the differences between CDM and 
a/m = 0.1 cm 2 /g are more subtle. 

To make a quantitative comparison between our simulations 
and the LoCuSS observations, we examine only those LoCuSS 
clusters that have ao and r core in a similar range as the fits to 
the five most massive halos in our simulations. This restricts the 
number of relevant LoCuSS cluster to five. In Fig. [TT| we show 
the ellipticity e distribution of the five most massive clusters in the 
CDM (black), a/m — 1 cm 2 /g (cyan), and a/m — 0.1 cm 2 /g 
(green hatched) simulations for lines of sight along the three prin- 
cipal axes of the moment-of-inertia tensors. The dark blue points 
with error bars show the central values and l-a uncertainties in 
e for the cluster halos in the lens-model fits for the five similar 
LoCuSS clusters. As expected, the ellipticities are highest for the 
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Figure 9. Top row: Surface density of a simulated galaxy cluster of mass M vir = 1.8 X 10 14 M Q in CDM (left), a/m = 0.1cm 2 /g (center) and 
<j/m = 1 cm 2 /g (right) cosmologies. The central regions in which the projected radius R < r m i n are masked as they are in the fits. These surface densities 
result from viewing the halo along the major axis of the moment-of-inertia tensor of all particles in the halo. Bottom row: Best-fit dPIE surface-density fits for 
the surface densities above. The ellipticities of the dPIE fits are: CDM, e = 0.28; a/m = 0.1 cm 2 /g, e = 0.25; a/m = 1 cm 2 /g, e = 0.29. 



intermediate-axis line of sight. In this instance, the ellipticities of 
the CDM and a/m — 0.1 cm 2 /g halos are more consistent with 
the LoCuSS sample than the a/m = 1 cm 2 /g halos. However, the 
lensing probability is highest for lines of sight closely aligned to 
the major axis jvan de Ven, Man delbaum & Keeton 2009; |Man-| 
delbaum, van de Ven & Keeton 2009), and in this case even the 
CDM halos appear slightly rounder than the LoCuSS clusters. The 
a/m = 1 cm 2 /g halos are definitely too round. 

Based on this initial set of LoCuSS halos, we believe it is 
safe to say that a/m — 0.1 cm 2 /g is at least as consistent with 
observations as CDM, but that there is significant tension with 
a/m = 1 cm 2 /g. There are several things that preclude us from 
making any statements stronger than this. First, we have a small 
sample of both simulated and observed galaxy-cluster halos. Sec- 
ond, we do not know the virial mass or alignment of the LoCuSS 
galaxy clusters, nor do we have a good handle on the selection func- 
tion of the survey. This means that we cannot make a direct compar- 
ison between the simulations and the observations. Third, we do not 
have any galaxy-cluster-mass halos with A/ vir > 2.2 x 1O 14 M 
in our simulations, and thus our ellipticity probability distributions 
for the lowest-tJo LoCuSS clusters are almost certainly biased, al- 
though it is not clear in which direction that bias goes. Based on 
the X-ray data, it appears likely that several of the LoCuSS clusters 
in our subsample are more massive than any simulated cluster even 
though crn and r core are similar to the simulated clusters ( |Richard| 
|etal.|2010) . Given that the dPIE fits are made based on the very cen- 



tral regions of the clusters, it is not unexpected that virial masses 
of the clusters can be very different even if the inner regions look 
similar. Finally, we do not include baryons in our simulations; it re- 
mains unclear how the presence of baryons alters the density profile 
and shape of galaxy clusters ( Scannap ieco et al.|2012| l. 

However, it is fair to say that this ensemble of observed galaxy 
clusters already places stronger constraints on the SIDM cross sec- 
tion than the other constraints we considered in this section, in par- 
ticular the constraint from MS 2137-23. While a/m — 1 cm 2 /g 
is easily allowed by our reanalysis of the MS 2137-23 constraint 
(modulo the uncertainty in the normalization of the convergence), 
this value of the SIDM cross section is in some tension with the 
LoCuSS cluster sample. A more quantitative limit, though, will 
only be possible with better theoretical predictions for the shapes of 
cluster-mass halos and a careful analysis of observed cluster selec- 
tion functions. In the future, we will simulate larger dark-matter ha- 
los and perform mock observations of them to find a better quantita- 
tive mapping between observations and SIDM cross-section limits. 

It may be difficult to probe cross sections as small as a/m — 
0.1 cm 2 /g based on cluster halo shapes alone, though. We get a 
sense of how hard it may be to use strong lensing to probe halo 
shapes, and hence small self-interaction cross sections, in Fig. [9] in 
which the shapes of the simulated CDM and a/m — 0.1 cm /g 
halos do not look very different on the scales to which strong 
lensing is sensitive. Moreover, in Rocha et al. (2012), we esti- 
mated that core sizes for cluster-mass halos should be < 20 kpc if 
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Figure 10. Same as Fig.[9lbut with the line-of-sight along the intermediate axis of the halo. The ellipticities here are: CDM, e = 0.59; a/m = 0.1 cm 2 /g, 
1 cm 2 /* 
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Figure 11. Ellipticity e (Eq. [8j of halos fit with dPIE profiles for three different projections. The solid black histograms show e values for the five most 
massive CDM halos in the 50 /i _1 Mpc simulation, the cyan histogram shows the same halos in the SIDMi simulation, and the green hatched histogram is for 
SIDMo.i. The dark blue points with uncertainties show the best-fit ellipticities and their 1-<t uncertainties from dPIE modeling of the five LoCuSS clusters 
with ctq and r core similar to those of the simulated clusters (Richard et al .|2010} . 



SIDM Halo Shapes 15 



a/m — 0.1 cm 2 /g. These sizes are similar to the effective radii of 
the brightest cluster galaxies (BCGs) at the centers of halos. There 
are several implications of this fact. First, it means that stellar kine- 
matics of the BGCs will be important for probing the dark-matter 
halo on scales for which scattering matters. Strong lensing is less 
sensitive to both the density profile and halo shapes at such small 
scales (see, e.g., Fig. 3 in |Newman et al.| (2011|). Second, it means 
that any inferences on the dark-matter halo on such scales depends 
on careful and accurate modeling of the BGCs in the data analysis. 
Third, we would have to model the behavior of SIDM in the pres- 
ence of a significant baryon-generated gravitational potential, and 
to explore the coevolution of dark-matter halos and BGCs. In par- 
ticular, simulations of isolated disky galaxies indicate that the pres- 
ence of baryons tends to make the dark-matter distribution more 
spherical, although it is not clear how much the dark-matter distri- 
bution changes if the central galaxy is elliptical instead (D ebattista| 
|et al.|20 08 1. We note that while these issues also have implications 
for SIDM constraints based on radial density profiles or central 
densities, they may be more serious for the shape-based constraints 
because the shapes of a/m — 0.1 cm 2 /g are already so similar to 
CDM halos for small radii even in the absence of baryons. 



5 CONCLUSIONS 

The takeaway message of this work is that mapping observations 
to constraints on the self-interaction cross section of dark matter 
is significantly more subtle than previously assumed, and as such, 
constraints based on halo shapes are, at present, one to two orders 
of magnitude weaker than previously claimed. 

There are three primary reasons that contribute to this con- 
clusion. First, the observational probes (gravitational lensing and 
X-ray surface brightness) of halo shapes are actually probes of 
some moment of the mass distribution. For lensing, the observa- 
tional probes are also sensitive to all material along the line of sight. 
While SIDM makes the three-dimensional density distribution sig- 
nificantly rounder within some inner radius r, the surface density 
will in general not be spherically symmetric at a projected radius 
R = r. The surface densities are affected by material well out- 
side the core set by scatterings where material is still quite triaxial 
(Fig. [3}. Previous constraints were made under the assumption that 
the observations tracked the three-dimensional halo shape for fixed 
projected radius. This is a less troublesome assumption for X-ray 
isophotes, since it is weighted by the square of the gas distribution, 
and hence sensitive to the central regions. The shapes measured 
should be related most closely to the shapes of the enclosed mass 
profile. So, to probe cross-sections as small as a/m = 0.1 cm 2 /g, 
one needs to get down to O(10 kpc) from the center of the halos. 
The contribution of stars in this region makes it difficult to robustly 
estimate the shape of the dark matter profile and it also makes it 
difficult to get a large ensemble of galaxies for this study. 

Second, there is a fair bit of scatter added by assembly history 
to the observed shapes and the scatter is large enough that it pre- 
cludes using a small number of objects to set constraints on SIDM 
cross sections. Finally, although we find that the three-dimensional 
shape of halos begins to become more spherical than CDM at radii 
where the local interaction rate is fairly low, T(r) ~ 0.1 Ho. there 
is a fair amount of triaxiality even when T(r) ~ Ho, a fact that 
was not appreciated in earlier studies (e.g., Miralda-Escud e]2002| 
|Feng, Kaplinghat & Yu|20"T0l > 

We find that the convergence map of MS 2137-23 (M v i r ~ 
1O 15 M0) allows a velocity-independent SIDM cross section of 



a/m = 1 cm 2 /g- The X-ray isophotes of NGC 720 (M vir ~ 
1O 13 M0) likely rule out a/m — 1 cm 2 /g, but are consistent with 
a/m — 0.1 cm 2 /g at radii where we can resolve shapes in our 
simulations. Based on a preliminary comparison to lensing mod- 
els of LoCuSS clusters, we conclude that a/m = 0.1 cm 2 /g is as 
consistent with observations as CDM but that a/m = 1 cm 2 /g is 
likely too large to be consistent with the observed shapes of those 
clusters. 

Cross sections in this range are very interesting. In Rocha et 
al. (2012), we show that a cross section in the neighborhood of 
a/m = 0.1 cm 2 /g could solve the "Too Big to Fail" problem 
for the Milky Way dwarf spheroidals (Boylan-Kolchin, Bullock & 



Kaplinghat 2012), the core-cusp problem in LSB galaxies (Kuzio 
|de Naray et al.|20T0)|de Blok|2010), as well as the shallow density 
profiles of the galaxy cluster s in|Sand et aT7| (2008 1; N ewman et al.| 
(2009), and New man et al.| (201 1 ) while not undershooting their 
central densities or overshooting the core sizes. Cross sections in 
this range are also consistent with other density-profile-based and 
subhalo-based constraints ( Yoshida et al. 2000b ; Gnedin & Ostriker 
[2001) . 

Since the current set of observations appear to be consistent 
with a SIDM cross section of a/m. = 0.1 cm 2 /g, there are two 
relevant questions for shape-based SIDM constraints. Will shape- 
based constraints be competitive with other types of SIDM con- 
straints? And what will it take to get down to a/m ~ 0.1 cm 2 /g 
with shapes? 

Upon closer inspection, our view is that constraints using ex- 
isting data could be pushed below a/m — 1 cm 2 /g, but it is 
not yet clear that we can get to a/m ~ 0.1 cm 2 /g. While X- 
ray isophotes of the elliptical galaxy NGC 720 are consistent with 
a/m = 0.1 cm 2 /g, there are some differences and a larger ensem- 
ble of elliptical galaxies may be able to test that. There are a number 
of other elliptical galaxies for which high-resolution X-ray data ex- 
ist (e.g., Humph rey et al.|2"0 06), but they lack the detailed shape 
measurements of NGC 720. So better constraints could result from 
X-ray shape analysis for these galaxies. For clusters, based on our 
quick pass through the LoCuSS results, lensing-based shape con- 
straints on SIDM could also extend well below a/m = 1 cm 2 /g 
if simulations are performed of a statistically significant number of 
massive galaxy clusters. However, in studies of both galaxies and 
clusters, it is likely that the measured densities in the inner regions 
would be a better way to test for signatures of self-interacting dark 
matter. 
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